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1.  Introduction 


Recently,  the  authors  developed  a  new  approach  for  the  study  of  generalized  stability 
of  families  of  real  (or  complex)  matrices  and  polynomials  using  so-called  guardian  and 
semiguardian  maps  (Saydy  et  al  (1990,  1988a,  1989)  and  Saydy  (1988b)). 

In  this  paper,  after  recalling  some  of  the  main  results,  we  present  a  closed-form  expres¬ 
sion  for  the  maximal  range  of  stability  of  one-parameter  families  of  matrices  or  polynomials 
relative  to  many  domains  of  practical  interest  of  the  complex  plane.  The  determination 
of  the  maximal  range  requires  finding  the  zeros  of  a  polynomial  which  depends  on  the 
family  and  domain  of  stability  under  consideration.  This  framework  is  then  applied  to  the 
study  of  three  physically  motivated  examples:  (1)  crane  stabilization  with  damping  and 
stability  margin  specifications,  (2)  satellite  attitude  control  and  (3)  control  of  a  digital 
tape  transport. 

2.  Guardian  Maps  and  Robust  Stability 

The  guardian  map  approach  was  introduced  in  Saydy  et  al  (1988a,  1988b,  1990)  as  a 
unifying  tool  for  the  study  of  generalized  stability  of  parametrized  families  of  matrices  or 
polynomials.  Some  of  the  basic  concepts  used  in  this  approach  are  reviewed  below.  We 
use  the  notation  T>  and  dT>  to  denote  the  closure  and  the  boundary  of  a  given  set  T>. 

2.1.  Guardian  Maps 

Basically,  guardian  maps  are  scalar  valued  maps  on  the  set  of  real  nxn  matrices1  that 
take  nonzero  values  on  the  set  of  “stable”  matrices  and  vanish  on  the  boundary  of  that  set. 
As  the  concept  of  guardian  maps  is  normally  applied  to  to  generalized  stability  problems 
(wherein  eigenvalues  are  confined  within  open  subsets  of  the  complex  plane  other  than  the 
open  left-half  plane),  it  is  useful  to  allow  the  set  of  stable  matrices  to  be  any  given  open 
subset  of  RTlXn.  Thus  we  have  the  following  definition. 

1  While  in  this  paper  we  restrict  ourselves  to  families  of  matrices,  all  results  extend 
readily  to  families  of  polynomials  as  well  (Saydy  et  al .,  1990). 
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Definition  1.  (Saydy  et  al,  1990)  Let  5  be  an  open  subset  of  HraXn  and  let  u  map  ]RnXn 
into  C.  We  say  that  u  guards  S  if  for  all  A  £  <5,  the  equivalence 

v(A)  —  0  -£=>•  A  £  dS  (1) 

holds.  The  map  v  is  said  to  be  polynomic  if  it  is  a  polynomial  function  of  the  entries  of  its 
argument. 

Of  special  interest  are  sets  of  the  form  5(0),  where 


5(0)  =  {ie  Tr"  :  <r(A)  C  0},  (2) 

O  is  an  open  subset  of  the  complex  plane  and  cr(A)  denotes  the  spectrum  of  A.  Such  sets 
5(0)  will  be  referred  to  as  ( generalized )  stability  sets. 

For  example,  the  map  u  :  A  ■— >  det  (A  ©  A),  where  ©  denotes  Kronecker  sum,  guards 

O 

the  set  ofnxn  Hurwitz  stable  real  matrices  5(C_).  This  follows  from  the  property  that 
the  spectrum  of  the  Kronecker  sum  of  two  square  matrices  A  and  B  consists  of  all  pairwise 
sums  of  eigenvalues  of  A  and  B  (Lancaster  and  Tismenetsky,  1985). 

2.2.  Robust  stability 

The  robust  stability  problem  for  parametrized  families  of  matrices  or  polynomials  may  be 
stated  as  follows.  Let  r  —  (r\ , . . . ,  r*)  £  U,  where  U  is  a  pathwise  connected  subset  of  IR*, 
and  let  A(r)  be  an  element  of  !RnXn  which  depends  continuously  on  the  parameter  vector 
r.  Given  an  open  subset  5  of  IFtnXn,  we  seek  basic  conditions  for  A(r)  to  lie  within  5 
for  all  values  of  r  in  U.  The  next  theorem  gives  a  basic  necessary  and  sufficient  condition 
for  this  problem  both  for  guarded  sets  5.  Typically,  5  is  a  stability  set  of  the  form  5(0) 
where  0  is  a  given  subset  of  the  complex  plane. 

Theorem  1.  (Saydy  et  al ,  1990)  Let  5(0)  be  guarded  by  the  map  v.  The  family 
{A(r)  :  r  £  U}  is  stable  relative  to  0  if  and  only  if 

(i)  it  is  nominally  stable,  i.e.,  A(r°)  £  5(0)  for  some  r°  £  U ,  and, 

(ii)  v(A{r))  0,  Vr  £  U. 
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In  the  case  of  polynomic  guardian  maps,  the  theorem  above  yields  computable  con¬ 
ditions  for  robust  stability  of  polynomially  parametrized  families  of  matrices.  In  Saydy  et 
al  (1990),  polynomic  guardian  maps  are  constructed  for  many  stability  sets  of  practical 
interest. 

2.3.  Maximal  interval  of  generalized  stability 

Let  Q,  be  an  open  subset  of  the  complex  plane  such  that  is  guarded  by  a  polynomic 
map  v  and  consider  the  one-parameter  family  of  matrices 

A(r)  =  A0  +  rAi  +  . . .  +  rmAm,  (3) 

where  A0,. . . ,  Am  are  given  real  square  matrices.  Let  r°  be  a  nominal  value  of  interest 
and  assume  that  A(r°)  is  stable  relative  to  Q. 

Denote 

u(r)  :=  u(A(r))  .  (4) 

We  seek  to  find  the  largest  open  interval  of  parameter  values  (rmj„,rmax)  containing  r°  for 
which 


A(r)  £  S(Q)  V  r  £  ( 

Dnin  ?  Dnax  )•  (5) 

Note  that  since  «S’(fi)  is  open  and  A(r)  is  continuous  in  r,  the  fact  that  A(r°)  £  <S(D) 
guarantees  the  existence  of  an  open  interval  containing  r°  in  which  the  family  above  is 
stable. 

The  theorem  below  easily  follows  from  Theorem  1. 

Theorem  2.  Let  <S(Q)  be  guarded  by  a  polynomic  map  u  and  let  A(r°)  be  stable  relative 
to  fi.  Let2 

Z~  =  {r  <  r°  :  u(r)  =  0}  U  {— oo}  (6) 

Z+  =  {r  >  r°  :  u(r)  =  0}  U  {+oo}.  (7) 

2  Note  that  i /(r°)  ^  0  since  A(r°)  £  S(il)  and  S( fl)  is  guarded  by  v. 
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Then 


v min  =  max  {z  :  z  €  Z  } 

(8) 

rmax  =  min  {z  :  z  G  Z+j 

(9) 

n 

Note  that  in  the  case  of  Hurwitz  stability  of  the  convex  hull  of  two  matrices,  it  was  shown 
in  Fu  and  Barmish  (1988),  based  on  a  result  of  Bialas  (1985),  that  the  maximum  range 
can  be  obtained  from  the  eigenvalues  of  a  suitably  constructed  matrix. 

3.  Examples  of  Application 


We  apply  the  results  of  the  previous  section  to  three  examples.  The  last  two  from  Franklin 
et  al  (1986)  were  studied  in  Bhattacharyya  (1986)  (see  also  Biernacki  et  al.  (1987))  where 
controllers  achieving  robust  stabilization  within  a  prescribed  set  in  the  the  parameter  space 
were  synthesized. 

Example  1:  Crane 

In  this  first  example  we  treat  the  case  of  a  crane,  a  simplified  model  of  which  consists 
of  two  masses  connected  by  a  light  inextensible  rod  (Fig.  1).  This  system  may  be  described 
in  state  space  form  by  (Hwang  and  Schmitendorf  (1984)) 


where 


x(t )  =  A(r)x(t)  +  b(r)u 


A(r ) 


(01  0  0\ 

0  0  r  0 

0  0  0  1  I  ’ 

0  0  — (1  +  r)  0/ 


6(r) 


and  r  :=  jfj  is  an  uncertain  parameter. 

Our  goal  is  to  synthetize  a  feedback  controller  u  =  kx  which  stabilizes  the  above 
uncertain  system  around  the  nominal  value  r°  =  1  with  full  information  on  the  maximal 
range  of  stability  as  well  as  guaranteed  damping  ratio  and  stability  margin.  To  this  end, 
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let  us  assume  that  it  is  desired  to  place  the  eigenvalues  of  the  above  uncertain  system  in 
the  domain  Cl  of  Fig.  2  specified  by  a  slope  a  =  2  and  a  margin  a  =  0.5  with  complete 
knowdelge  of  the  maximal  parameter  range  for  which  these  eigenvalues  remain  in  Cl.3 


Fig.  1.  Crane  of  Example  1 


Fig.  2.  Stability  domain  for  Example  1 


Note  that  the  controllability  of  the  nominal  system  makes  it  possible  to  prespecify  Cl. 
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Accordingly,  let  A  =  {— 1  ±  0.5j,  —  1  ±  j}  C  ft  be  desired  closed- loop  eigenvalues  for 
the  nominal  system.  Then  one  easily  finds  that  the  required  feedback  is  [1.25  3.25  — 
1.375  1.25]  :=  ko.  To  determine  (rmin5  rmax),  the  largest  range  of  stability  relative  to  ft, 
we  may  apply  Theorem  2  to  Aci(r)  :=  A(r)  —  bko ,  provided  we  can  endow  <S(ft)  with  a 
polynomic  guardian  map.  It  turns  out  that  the  mapping  given  by  z/(A)  =  iq(A)i/2(A), 
where 

ui(A)  =  det  |  ^ —A  ■  A  -  A2  •  j| 

and 

^(A)  —  det  ((A  +  crl)  •  I)  det  (A  +  al)  , 
is  one  such  map  (Saydy  et  a J,  (1990)).  We  obtain  the  polynomials 

vi(Ad(r))  =  1.933  104(r6  +  4.92757  r5  +  9.04694r4  +  6.94471  r3  +  0.828279r2 
-  1.62776  r-  0.630638) 

+  80r  -  25)  ^  (9r  +  11) 

which  vanish  at  the  real  values 

{-1.0,  -1.4399487,  -0.9490848,  0.4614538252), 

{-9.927086,  -0.670474,  0.20867,  -1.22222), 

respectively.  Recalling  that  r°  =  1,  we  conclude  by  virtue  of  Theorem  2  that  rmjn  = 
0.4614538252,  rmax  =  +oo,  i.e.  that  the  closed  loop  uncertain  system  is  asymptotically 
stable  for  all  values  r  G  (rm;n,  oo)  and  unstable  at  r  =  rm-m.  (At  this  value  of  r,  the  closed 
loop  eigenvalues  are  -0.76734499  ±j  0.4117427  G  ft,  and  -0.694108944  ±j  1.388217872  G 
5ft.) 

Example  2:  Satellite  attitude  control 

A  satellite  system  is  modeled  as  two  masses  connected  by  a  spring  with  torque  constant  k 
and  viscous  damping  constant  d  (Fig.  3).  The  equations  of  motion  are  given  by 

J\9\  +  d{6\  —  62)  +  k(9\  —  9 2 )  =  Tc 

J2&2  A  d(9 2  —  9\ )  -f-  k{92  —  9\ )  =  0 
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Model 


Fig.  3.  Sketch  of  a  satellite 


where  J\  and  J2  are  inertias  and  Tc  is  the  control  torque.  Choosing  the  state  vector  to  be 
x  —  and  letting  u ,  r\,  r2  denote  denote  Tc ,  k  and  d  respectively,  we  obtain 

the  state  equations 

x  =  A(ri ,  r2)x  +  bu 


y  =  cx 


with 

0  10  0 

_  7*1  7*2  ri  7*2 

J 2  J 2  J2  J? 

0  0  0  1 

_£l_  _£2.  £L  £2. 

Jl  Jl  Ji  J\ 


c  =  (0,0, 1,0). 


A(n,r2)  = 
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Physical  analysis  leads  to  the  conclusion  that  the  parameters  n  and  r2  vary  as  a  result  of 
temperature  fluctuations  within  the  range  (Franklin  et  al. ,  1986): 


0.09  <  n  <  0.4, 


0.04 


<  r2  <  0.2 


The  matrix  A(ri,r2)  is  singular  for  every  value  of  the  parameters  ri  and  r2  and  is  hence 
Hurwitz  unstable.  Letting  u  =  —ky  for  some  scalar  k ,  we  obtain  a  closed  loop  matrix 
A *(rj,r2)  :=  A(ri,r2)  -  k  b  c. 

O 

As  we  have  seen,  the  set  of  Hurwitz  stable  matrices  <S(C_)  is  guarded  by  the  polynomic 
map  v  :  A  det(A  ®  A).  Rather  than  using  this  map  (it  involves  in  our  case  the 
computation  of  a  16  x  16  determinant),  we  use  the  alternate  (polynomic)  guardian  map 
(see  Propositions  7  and  8  in  Saydy  et  al  (1990)) 


y;  Ah  ui(A)v'2(A) 


with 

vx(A)  =  det(A  •  I) 
i/2(A)  =  det(A) 

where  the  matrix  A  •  I  denotes  the  bialternate  product  of  A  and  I  and  has  dimension  6. 4 
With  J\  =  J2  =  1  we  obtain5 


*/(^*(n,r2))  =  (g4^3r2)  (n) 


^  „2 
r.r2. 


We  can  therefore  conclude  using  Theorem  1  that  if  k  is  any  controller  which  stabilizes 
the  system  for  some  nominal  values  r°  >  0,  r2  >  0,  then  it  also  stabilizes  the  satellite 

4  For  an  n  x  n  matrix  A ,  the  matrices  A  ©  A  and  A  •  In  are  n2  x  n2  and  x  H— ^12 

respectively. 

5  Computer  algebra  codes  were  used  to  carry  out  the  computations  symbolically  in  both 
examples. 
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system  globally  in  the  open  first  quadrant,  i.e.  for  all  strictly  positive  r\  and  r2,  but 
stability  is  lost  if  n  or  r2  vanishes. 

Example  3:  Control  of  a  digital  tape  transport  system 

The  objective  in  this  last  example  is  to  control  the  speed  and  the  tension  of  the  tape  at 
the  read/write  head.  A  model  of  such  a  system  is  (Bhattacharyya,  1987) 

x  =  A(r\,r2)x  +  b  u 
y  =  c  (rur2)x 

where  A(r\ ,  r2)  is  a  5  x  5  companion  matrix  with  last  row  given  by 
A5i(r1,r2)  -  0 
^52(ri,r2)  =  -0.1045  10~4n 
A53(ri,r2)  =  -(0.35  lO^n  +  0.01045  r2  +  0.2725) 

A54(rur2)  =  -(0.25  10~4n  +  0.035  r2  +  1.5225) 

A55(rur2)  =  -(0.025  r2+2.25) 

and  where 

b  —  (0, 0, 0, 0,  l)r,  c(n,r2)  =  (0.3  10-4n,  0.03  r2,  0,  0,  0). 

It  is  assumed  that  the  parameters  r\  and  r2  are  subject  to  perturbations  around  their 
nominal  values  of  r9  =  4.0  104  and  r°  =  20.  Let  Ak(ri,r2)  again  denote  the  closed  loop 
matrix  A(ri,r2)  —  bkc(ri,r2).  It  is  shown  in  (Bhattacharyya  (1987),  p.  71)  that  the  gain 
k  =  0.00827  stabilizes  the  system  robustly  with  stability  margin  9.99.  We  propose  to 
analyze  the  performance  of  this  controller  using  the  approach  presented  here.  Using  the 
guardian  map  of  the  previous  example  we  obtain  for  k  —  0.0827  that 

V!  (Ak(rur2))  =  1.5258  10“23[1.2056  1012  r\  +  (5.0782  1  09  n  +  1.7002  1014)r3 

+  (1.7008  107  r\  +  3.9995  10n  rx  +  7.8397  1015)r^ 

+  (1.214  104  r3  +  1.4755  109  r\  +  1.6002  1012n  +  1.3643  1017)r2 

+  (4.1306  105  r3  +  3.1718  104  r2  -  1.8761  1014  n)] 

v2{Ak(ri,r2))  =  -2.48  10-6  rx. 
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Let  us  freeze  one  of  the  parameter  at  its  nominal  value,  say,  r\  =  rj  =  4.0  104. 
Then  one  can  show  using  Sturm  sequences  for  example  that  the  univariate  polynomial 
u\{Ak{r\.,  r2))u2(Ak(r^  >  ^2))  has  no  zeros  in  [0,  +00),  implying  by  virtue  of  Theorem  1  that 
the  family  r2)  is  Hurwitz  stable  for  all  values  r2  £  [0, 00).  Repeating  the  same 

computation  in  the  direction  r2  =  r%  =  20  yields  the  similar  conclusion  that  the  fam¬ 
ily  Ak(ri,r2 )  is  Hurwitz  stable  for  all  n  £  (0, +00).  This  suggests  that  the  controller 
k  =  0.0827  possibly  achieves  stability  for  all  positive  values  of  the  parameters  r\,  r 2. 

To  investigate  this  question  we  use  the  following  result: 

Fact.  (Saydy  et  al. ,  1990)  Let  f2  be  a  subset  of  the  complex  plane  such  that  <S(f2)  is  guarded 
by  a  given  real  polynomic  map  v.  Let  A  {4.(ri,r2)  :  (r\,r2)  £  [oi,/?i]  x  [0:2, #2]}  be 
a  nominally  stable  family  of  real  matrices;  e.g.,  .4(0:1,02)  £  <S(fi).  Then  the  family  A 
is  stable  relative  to  Q  if  and  only  if  U*T  =  0  and  the  univariate  polynomials  pai ,  and 
pr i)  r\  £  U^r,  have  no  zeros  in  [02,  /32\.  Here,  for  each  ri,  pri  denotes  the  univariate 
polynomial  v(A(r\ ,  •)), 


^cr  :=  {ri  €  [01  ,/?i]  :  det  B{pn  ,p'ri)  =  0}  , 

^cr  :=  {n  e  [a i,/?i]  :  Pr1(o2)Pn(/3 2)  =  0}  . 

where  JE?(p,  q)  denotes  the  Bezoutian  of  the  polynomials  p  and  q  and  prime  denotes  deriva¬ 
tive. 

□ 

Note:  the  following  remarks  regarding  the  application  of  the  fact  above  are  in  order: 

(i)  the  interval  within  which  rx  lies  may  be  any  interval  (i.e.  not  necessarily  closed); 

(ii)  the  factor  pn(/? 2)  may  be  omitted  in  U%r  if  /32  =  +00; 

(iii)  it  is  necessary  to  verify  that  the  polynomial  p0l  has  no  zeros  in  [a2,/32]  only  if  U*t 
turns  out  to  be  empty. 

Clearly  v2(Ak(rx ,  r2))  does  not  vanish  in  the  first  open  quadrant  and  we  can  thus  rede¬ 
fine  pri  as  vi(Ak(ri,r2)).  It  can  be  verified  using  Sturm  sequences  that  the  polynomial 
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det B(pn,p'n),  given  (modulo  a  constant  factor)  by 

1.7383  r\2  -  2.5952  104  rj1  -  4.3427  109  r\°  +  2.0134  1014 

-  1.3604  1018  r8  -  1.403  1023  r\  +5.2263  1  027  -  8.5599  1031 

+  6.5237  1035  r\  -  2.9793  1038  rf  -  2.4777  1043  rj  +  4.896  1046  n  +  5.1342  105° 

has  only  one  zero  in  (0,  -f-oo),  namely  14005.4616.  With  the  notation  of  the  fact  above  and 
(ri,^)  €  (0, +oo)  x  [0,+oo)  we  thus  have 

UlT  —  {ri  €  (0,oo)  :  det  B(pri,p'ri)  =  0}  =  {14005.4615} 

U2r  =  {ra  e  (0,  oo)  :  pn (0)  =  0}  =  {5518.349}. 

It  follows  that  the  closed-loop  matrix  Ak(ri,r2)  is  not  Hurwitz  stable  for  all  (ri,r2)  in 
(0,oo)  x  [0,  oo)  since  £/2r  ^  0.  In  fact  one  can  easily  check  that  for  ri  =  14005.4615,  the 
polynomial  pri  has  no  zeros  in  [0,  oo).  We  therefore  conclude  by  virtue  of  the  fact  above 
that  the  Hurwitz  stability  of  the  closed-loop  system  is  guaranteed  in  the  smaller  rectangle 

(n^)  €  (5518.349, +oo)  x  [0, +oo). 
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